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Abstract 

The fluctuations and the distribution of the conductance peak spac- 
ings of a quantum dot in the Coulomb-blockade regime are studied and 
compared with the predictions of random matrix theory (RMT). The ex- 
perimental data were obtained in transport measurements performed on 
a semiconductor quantum dot fabricated in a GaAs-AlGaAs heterostruc- 
ture. It is found that the fluctuations in the peak spacings are considerably 
larger than the mean level spacing in the quantum dot. The distribution 
of the spacings appears Gaussian both for zero and for non-zero magnetic 
field and deviates strongly from the RMT-predictions. 

PACS numbers: 73.20.Dx,73.23.Hk,05.45.+b 

Advanced nanofabrication techniques have made it possible to confine small 
numbers of electrons electrostatically within the two-dimensional electron gas 
(2DEG) of a semiconductor heterostructure pi ||. Both the electric charge and 
energy of such "quantum dots" are quantised and hence such structures are 
sometimes referred to as "artificial atoms" || ||. In transport measurements 
the charging of these electron islands with single electrons leads to the obser- 
vation of periodic conductance oscillations in the Coulomb-blockade regime jl| . 
These reflect the electrostatic coupling of the quantum dot to its environment 
and, additionally, they contain information about the eigenenergies and eigen- 
functions of the electrons in the dot. Due to irregularities in the electrostatic 
confinement potential and electron-electron interactions, the corresponding clas- 
sical motion of the electrons in the quantum dot can be expected to be chaotic 
(nonintegrable) || [(], . Consequently, recent experiments have considered the 
peak height distribution || ||] , parametric conductance correlations || and level 
statistics [[To] of a quantum dot in the Coulomb-blockade regime to test the con- 
cepts developed for the quantum mechanical description of classically chaotic 
systems ( "quantum chaos" JLl| [12[ ) . In particular random matrix theory (RMT) 
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[13] has proven to be a very successful description of the statistical properties of 
spectra of many irregular systems. Therefore, it is a very interesting question, 
how RMT applies to the transport properties of quantum dots. In this paper 
we investigate the fluctuations of the peak spacings of the conductance peaks of 
a quantum dot obtained in transport measurements with and without applied 
magnetic field. The spacing distributions are calculated and compared with the 
predictions of RMT. 

For sufficiently low temperatures T and small dot capacitances C, quantum 
dots isolated from the reservoirs of the 2DEG via tunnel barriers can exhibit 
Coulomb-blockade phenomena. When e 2 /C ^> fc^T, transport through the 
quantum dot is blocked. A finite conductance only occurs when the total en- 
ergy of the quantum dot with N electrons is degenerate with the energy of the 
dot occupied by N + 1 electrons. This is the case when 



where F(N) denotes the free energy of the quantum dot with N electrons and 
/i is the chemical potential of the leads. Then a single electron can tunnel from 
a reservoir into the dot jjj. This can be achieved by tuning the dot's potential 
with a centre gate. A sweep in the centre gate voltage V g results in the well- 
known conductance oscillations in the Coulomb-blockade regime. From Eq.(Q) 
the difference AV g between gate voltages at which two adjacent peaks occur 
can be related to the thermodynamic quantity d[i/dN, which has the meaning 
of an inverse compressibility [0 . Within the capacitive charging model [Q the 
electrons are assumed to occupy single particle states of energies e$ and the 
Coulomb interactions are described by a classical electrostatic term U(N). The 
dot's energy is then F(N) « £f + U(N) and the difference AV g is given by 



Here e denotes the electronic charge, C the total capacitance of the dot and 
Aejv = ejv+i — ejv the level spacing. The conversion factor a = C g /C, where C g 
is the dot-to-gate capacitance, translates between the energy and the voltage 
scale of the conductance oscillations. Thus, in principle, one should be able to 
extract the energy level spacings AeAr from the so-called "addition spectrum" 
obtained in Coulomb-blockade measurements. 

From the addition spectrum, one can calculate the nearest neighbour spacing 
(NNS) distribution P(S), which can be compared to the predictions of RMT. 
P(S) is the distribution of the spacings between adjacent levels of an energy 
spectrum, where the spacings S are normalised to a mean value of unity. The 
results for P(S) within the framework of RMT are very well approximated by 
the Wigner surmise, which is E2] 



F(N + 1) - F(N) = n, 



(1) 



eaAV g N = e 2 /C + Ae N . 



(2) 



P(S) 



P(S) 
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(GOE) 
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for time-reversal invariant systems and for systems with broken time-reversal 
invariance, e.g. in the presence of a magnetic field. The first distribution cor- 
responds to the energy spectrum of Hamiltonians drawn from the Gaussian or- 
thogonal ensemble (GOE) of random matrices, while the second is obtained for 
the Gaussian unitary ensemble (GUE). The fluctuations SS = ((5 2 ) - {S) 2 ) 1 ' 2 
are thus expected to be 0.52(5) and 0.42(5) for GOE and GUE, respectively. 
The quantum dot on which our measurements were performed was defined by 
electron-beam lithography in the 2DEG of a GaAs-Alo.32Gao.6sAs heterostruc- 
ture. The mobility and the sheet density of the 2DEG are 120m 2 /Vs and 
3.6 x 10 15 m~ 2 , respectively. The application of negative gate voltages to the 
surface structure defines an island which is isolated from the left and right 
reservoirs via tunnel barriers (see inset of Fig. [I]). The radius of the island is 
estimated to be r < 150nm. The Coulomb-blockade measurements were per- 
formed in a dilution refrigerator with a base temperature of 25mK. Electron 
transport through the dot was studied by applying a small bias voltage (4.3/^V 
AC) between the left and right reservoirs and measuring the current using stan- 
dard lock- in techniques (for further details see Ref. Bj). From the onset of 
the conductance oscillations at V g = — 200mV roughly 170 peaks are observed 
(Fig. [j]). During the gate sweep from V g = — 200mV to — lOOOmV the number 
of electrons in the quantum dot thus varies from N w 250 to N « 80. In the 
following, the only relevant energy scale is the mean energy level spacing A. It 
should be roughly Ep/N, where Ep is the Fermi energy. From the sheet density 
one obtains Ep « 12.9meV and therefore A ss 50/ieV. The thermal energy fc^T 
is about one order of magnitude smaller. 

To calculate the NNS distribution from the conductance oscillations, first the 
gate voltage differences AVg between adjacent peaks are extracted from the 
data. The mean value of AV g increases linearly with decreasing voltage (se e 
Fig. ||(a)) reflecting an inverse linear change in the dot-to-gate capacitance pH . 
Identifying (AVg) with the classical charging voltage e/C g , from Eq.(||) the 
energy spacings are obtained as 

Ae = ea(AV g - (AV g )). (5) 

The conversion factor a is also a function of the gate voltage. This can be consid- 
ered by using the same linear fit as above, i.e. a = C g /C = e(e+C res f (AVg)) -1 , 
where the capacitance C res t = C — C g is assumed to be constant. However, the 
actual choice of a is not a crucial parameter in the calculation, as tests with 
different constant values for a have shown. From Eq.(^J) the Ae are obtained as 
fluctuations around a mean value of zero. To remove the unphysical negative 
values for Ae the whole data are shifted by a constant value (cf. Fig. ||(b)). It 
turns out, that the fluctuations around the mean value are considerably larger 
than the mean level spacing A estimated above. This may be regarded as an 
indication that the calculated Ae are not the "real" addition energies and that 
the influence of the electron-electron interactions both within the dot and its 
environment play a significant role |i0|. 
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From the energy spacings one can construct an artificial one-particle energy 
spectrum via Ei — Y?n=i ^n- To unfold the data to a mean level spac- 
ing of unity a polynomial fit is made to the spectral step function N(E) = 
9(E — Ei). The renormalized energies are obtained by the standard unfold- 
ing mapping Ei i— > (N(Ei)) O, |l2). From these, the energy level spacings are 
calculated and can be directly compared to the predictions of RMT. 
In Fig. ||(a) the resulting NNS distribution for zero magnetic field is displayed. 
It obviously does not agree with the Wigner surmise (Eq.(||)). Instead, it is 
much better described by a Gaussian centered at S = 1, as illustrated. 
In the presence of a magnetic field, time-reversal invariance breaks down. In this 
case the appropriate ensemble of random matrix theory to describe energy level 
fluctuations is the Gaussian unitary ensemble. However, as in the B = OT-case, 
the experimentally obtained spacing distributions look Gaussian rather than 
GUE-like. In Fig. |(b) the distributions for zero, for low [B = 0.1T,B = 0.5T) 
and high (B — 4T) magnetic fields are displayed. The distributions are de- 
rived with typically 150 — 170 data points. From a statistical point of view this 
number is rather small. Nonetheless, these are large numbers when compared 
to previous Coulomb-blockade experiments [[To) . It can be seen that the distri- 
butions narrow with increasing magnetic field as may be expected due to the 
Landau quantisation jl6], |l7| . 

The largeness of the fluctuations indicates that the capacitive term e 2 /C in 
Eq.(|2|) undergoes even larger fluctuations than the energy levels themselves. 
Thus the Ae obtained above do not display the energy level spectrum of the 
quantum dot. However, this would not mean a failure of RMT, but an insuf- 
ficiency of the capacitive charging model. Eq.(Q) obviously cannot be used to 
get access to the bare energy level spacings of the quantum dot, when a larger 
range of gate voltages is considered. 

In a recent publication Sivan et al. fllpf argued that electron-electron interac- 
tions in the dot were responsible for the failure of RMT to describe the con- 
ductance peak spacing distribution. Their experiments and calculations lead 
to a Gaussian P{S) centered at S = 1 which is similar to our results. In 
terms of the charging energy the fluctuations obtained in our experiment are 
S(Ae) ~ 0.07 — 0.11e 2 /C, which is slightly smaller than in the work by Sivan 
et al. Their numerical calculations suggest that fluctuations in the quantity 
AV/(AV) converge to a "universal" value between 0.1 and 0.2 for strong elec- 
tronic interactions. Calculating this quantity from our data we arrive at 0.10, 
which is consistent with their finding. However, the influence of the capacitive 
coupling to the reservoirs has not been considered in their publication, which 
may also have a considerable influence on the fluctuation properties of the peak 
spacings. 

Finally, it has to be considered that RMT was initially developed to handle the 
statistical properties of excitation spectra of complex systems. The addition 
spectrum as obtained in Coulomb-blockade measurements, however, consists of 
the many-particle ground state energies of the quantum dot rather than excita- 
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tion energies. The comparison with RMT has been made under the assumption 
that the addition spectrum be equivalent to a single particle spectrum. Indeed, 
the excitation spectrum of the model used in jl(|, obeys RMT. Recently, this 
could also be shown for the excitation spectrum of the two-dimensional Hubbard 
model |l8). But it is not clear whether the results of RMT can be applied to 
ground state energy statistics. In this respect it is interesting to notice, that the 
peak height distribution for the conductance oscillations seems to be in accor- 
dance with RMT H P , whereas the parametric conductance correlations |)|, [j~9| 
quantitatively do not agree with RMT. 

In conclusion, we have investigated the statistics of conductance peak spacings 
obtained in Coulomb-blockade experiments with zero and non-zero magnetic 
field. In all cases the results do not agree with the predictions of random ma- 
trix theory. Instead, the nearest neighbour spacings appear to be Gaussian 
distributed around their mean value. It seems to be difficult to extract the 
bare energy levels when using a simple capacitive charging model. Therefore, 
our results include fluctuations in the electrostatic coupling with the environ- 
ment which are larger than the fluctuations in the quantum dot's energy level 
spectrum itself. Further theoretical and experimental work are required to un- 
derstand this central phenomenon in mesoscopic physics. 

We gratefully acknowledge stimulating discussions with R. Bcrkovits, H. Bruus, 
W. Hausler, J. P. Kotthaus, C. Stafford and S. Ulloa and financial support from 
the Deutsche Forschungsgemeinschaft. One of us (F.S.) also acknowledges sup- 
port from the Institute for Scientific Interchange Foundation. 

* present address: ETH Hoenggerberg HPF C13, CH-8093 Zurich 

References 

[1] H. van Houten, C. W. J. Beenakker, and A. A. M. Staring in Single 
Charge Tunneling, edited by H. Grabert, J. M. Martinis, and M. H. Devoret 
(Plenum, New York, 1991). 

[2] U. Meirav and E. B. Foxman, Semicond. Sci. Technol. 10 255, (1995) and 
references therein. 

[3] M. A. Kastner, Physics Today 46, 24 (1993). 

[4] R. C. Ashoori, Nature 379, 413 (1996). 

[5] R. A. Jalabert, A. D. Stone, and Y. Alhassid, Phys. Rev. Lett. 68, 3468 
(1992). 

[6] A. D. Stone and H. Bruus, Physica B 189, 43 (1993). 



5 



[7] H. Bruus and A. D. Stone, Phys. Rev. B 50, 18275 (1994). 

[8] A. M. Chang, H. U. Baranger, L. N. Pfciffcr, K. W. West, and T. Y. Chang, 
Phys. Rev. Lett. 76, 1695 (1995). 

[9] J. A. Folk, S. R. Patel, S. F. Godijn, A. G. Huibers, S. M. Croncnwctt, C. 
M. Marcus, K. Chapman, and A. C. Gossard, Phys. Rev. Lett. 76, 1699 
(1996). 

[10] U. Sivan, R. Bcrkovits, Y. Aloni, O. Prus, A. Aucrbach, and G. Ben- 
Yoseph, Phys. Rev. Lett. 77, 1123 (1996). 

[11] M. C. Gutzwiller, Chaos in Classical and Quantum Mechanics (Springer- 
Verlag, New York, 1990). 

[12] F. Haake, Quantum Signatures of Chaos ( Springer- Verlag, Berlin, 1990). 

[13] M. L. Mehta, Random Matrices (Academic Press, London, 1991), 2nd ed. 

[14] T. Hcinzcl, S. Manus, D. A. Wharam, J. P. Kotthaus, G. Bohm, W. Klein, 
G. Trankle, and G. Weimann, Europhys. Lett. 26, 689 (1994). 

[15] D. B. Chklovskii, B. I. Shklovskii, and L. I. Glazman, Phys. Rev. B 46, 
4026 (1992). 

[16] P. L. McEuen, E. B. Foxman, J. Kinaret, U. Meirav, M. A. Kastner, N. S. 
Wingreen, and S. J. Wind, Phys. Rev. B 45, 11419 (1992). 

[17] T. Hcinzcl, A. T. Johnson, D. A. Wharam, J. P. Kotthaus, G. Bohm, W. 
Klein, G. Trankle, and G. Weimann, Phys. Rev. B 52, 16638 (1995). 

[18] H.Bruus and J.-C. Angles d'Auriac, Europhys. Lett. 35 (5), 321 (1996). 

[19] H. Bruus, C. H. Lcwcnkopf, and E. R. Mucciolo, Phys. Rev. B 53, 9968 
(1996). 



6 



Figure 1: Conductance oscillations of a quantum dot in the Coulomb-blockade 
regime at zero magnetic field as a function of the centre-gate voltage V g . Roughly 
170 peaks are observed between V g = — 200mV and V g = — lOOOmV. The inset 
shows a schematic of the quantum dot. The shaded area denotes the 2DEG 
and the black areas indicate the gates with which the dot is defined. The lower 
middle gate is the centre-gate. 



Figure 2: (a) The peak spacings (AV^) extracted from Fig. 1 (indicated by dots) 
and a linear fit to them, (b) Energy spacings (Ae) calculated from the peak 
spacings and shifted to positive values. The straight line indicates their mean 
value and the broken lines indicate the largeness of the fluctuations expected 
from RMT. 



Figure 3: (a) NNS histogram calculated from the energy spacings of Fig. 2(b) 
after unfolding them to a mean value of unity. The full line denotes the GOE- 
prediction of RMT for P(S) and the broken line is a Gaussian fit centered at 
S = 1. (b) NNS distributions for different magnetic field strengths. The full and 
the broken lines denote the RMT-predictions for P(S) obtained in the Gaussian 
orthogonal and unitary ensemble, respectively. For clarity, lines have been used 
to diplay the distributions instead of histograms. 
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